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The method of hyperspherical "surface" functions has been apphed to the calcu- 
lation of eigenvalues and eigenfunctions of muonic molecular ions ^'"^Hedfi. Binding 
energies and nonradiative decay rates for the states of the total angular momentum 
L = 0, 1, 2 have been obtained. 
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I. INTRODUCTION 



The reasons for investigation of the charge-nonsymmetric muonic molecules like HeHjj, 
are as follows. First, a direct charge-exchange reaction from the ground-state muonic hy- 
drogen atom to helium nuclei is suppressed, the transfer proceeds through the formation of 
the molecule in the intermediate state. Hence, the kinetics of muons in media is defined to a 
large extent by the probability of this process. Indeed, the role of the formation of a muonic 
molecule in a charge-exchange reaction was confirmed in a number of experiments |[T|, 0, 

§, i- 

The measurement of the yield of 7-rays due to the decay of the Hedjji molecules p, 
IP revives interest in the investigation of this system. The experiment gives evidence of an 
additional nonradiative decay channel. This possibility was discussed in papers 0, [§, 

i- 

Next, eigenenergies of both usual molecules of media and muonic molecules, as it follows 
from the results of [|l^], [jlT|, are comparable. For this reason one can expect an active 



interaction of muonic molecules with media. 

As soon as the charge-nonsymmetric molecules are produced, the possibility of nuclear 
transitions arises. The investigation of a nuclear reaction at typical mesomolecular ener- 
gies has a fundamental importance due to absence of any experimental data on the strong 
interaction of charged particles in this energy range. 

Qualitatively properties of the HeHjji system are defined as follows. Coulomb interaction 
is not able to bind the systems under consideration due to the repulsion in the He^ + H 
channel. Only a 3— body resonant state can be formed. States like that are supported by 
the attractive polarization potential in the + He channel and therefore are clustered. 

The goal of this paper is to perform systematical calculations of energy levels and non- 
radiative decay rates of the ^''^Hedfi systems for all possible values of the total angular 
momentum. Only the transition to the channel with the Hefi atom in the ground state has 
been considered. Transitions from the molecular states with L = 1, 2 to the channels with 
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the He/j, atom in the 2s, 2p states will be suppressed due to exponentially small overlapping 
of the initial and final state wave-functions. 

The treatment of this problem met essential difficulties due to necessity to describe the 
Coulomb three-body system above the two-body threshold. Some approaches to describing 
these systems have been applied in [|IT|, 0, [g], 0. 

The approach, using the hyperspherical "surface" functions method ||12|, ||13| has been 
applied in this paper. The following advantages of this method in treating the posed prob- 
lems can be mentioned. The method operates with a discrete set of coupled one-dimensional 
differential equations. Physical boundary conditions for their solution can be easily formu- 
lated. Moreover, coupling of channels turns out to be rather small in our case and allows 
one to use the decoupled one-level approximation. It is worthwhile to mention the analogous 



calculation of the LiHjj and BeH/j molecules ||14|| . 

The article is organized in the following way. The description of the method will be 
given in the next section, section 3 contains numerical results, section 4 - discussion and 
conclusion. 



II. METHOD 



The Hamiltonian of three charged particles in the Jacobi variables is: 

if = -Ax,-Ay, + f ^, (1) 



where 



s=l 



m{mk + nij] 



Xi = J— : rirfc - r^j. 



_ / mi{mk + mj) ^ rrijTj +mkVk ^ 
y myrrii -|- nik + nij) rrij + rrik 



(2) 



= 2Z,ZJ (3) 
y m{mk + rrij) 
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Tj, rrii, Zi - coordinate, mass and charge of the i-th. particle. (me^), me^/ (2^^) have been 
used as length and energy units. Here m is arbitrary mass and was taken equal to the 
muonic mass. For definiteness muon, hydrogen nucleus and nucleus of the charge Z have 
been enumerated as particles with number 1, 2 and 3. 

Coordinates p, ai, 6'ihave been introduced by the relations 

ai 

Xi = pcos—, 

. OLi 

y, = psm—, 

„ (x^y.) W 

COS&i = , 

< ttj, 9i < TT. 

Below the notation Q will be used for an arbitrary pair a,, 6i. 

Since the systems have two heavy and one light particles, it is reasonable to assume that 
the main part of the total angular momentum is carried by the pair of heavy particles. This 
is reason why the following form of the solution of the Schrodinger equation has been used: 

*LM(x,y) = YLMi^,)^Lip,^). (5) 

Under these assumptions the Schrodinger equation for the states with total angular 
momentum L takes the form: 



where 



^ I d ^^^^^ d 1 d ^^^^ d L{L + 1) 

^ sirfiai dai ^ dat sinOidOi ^ dOi Acos'^ai 



Following [T^ and [^] "surface" functions p) can be introduced as finite solutions 

of the equation: 

[Ac-^E^ + A„(p)]^„(fi;p) = 0. (8) 

^ s=l 

Expanding the solution of the equation (j^) onto the set of the hyperspherical "surface" 
functions: 



one immediately comes to the system of one-dimensional equations 
cf 15 



E[Q™(P)^ + J-pQniip) - Pmip)]Unip) = 0, (lO) 



where 



d 

Qni{p) = M-Q^^i), (11) 

d d 

Pni{p) = ^^'^"l^^*)' (12) 
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^n(p) = —Kip)- (13) 
P 

( ■ I ■ ) means the integration on the hypersphere over dQ = sin^aidaidcosOi. 

One of the most complicated problems of this approach is the computation of Qni{p) and 
Pni{p), defined in (|l^) and (0). By this reason, the following exact expressions have been 
used: 

Qni = - K)~^{^n\^—\'fi), (14) 

Pm = ~{Q')m. (15) 

The form (P^), ( |T3| ) allow one to avoid the calculation of the derivatives of the "surface" func- 
tions on the parameter p and use only already known matrix elements Vni{p) and eigenvalues 
Aj(p) of equation (H). 

The variational approach has been applied to solve equation (|). The "surface" functions 
have been chosen as a linear combination of trial functions from the following set: 

(j)i"^{a^)Pi{cos9^), a = 2,3, 

szn^a^CltUPiicose^), (16) 
n> 0, n> l> 0, 



where 

(a) = ILii^-^pcos^), ^^^^ 

In equations (|16D and (|l^) Pi(x), C™(x) are the Legendre, Laguerre and Gegenbauer 

polynomials. The set of trial functions has been chosen in the form (p!6D in order to describe 
properly the three-body wave-function at both large and small interparticle distances. The 
first line of (p^ ) will describe the system separated into two clusters. In this case, one of 
the clusters is a hydrogen-like atom and hydrogen-like functions (p^Tf ) will be proper trial 
functions. The second line of (|16D will describe the configuration with all three particles 
close to each other. In this case, the centrifugal term in (P) dominates and eigenfunctions 
of the operator are used. The set of trial functions (|l6Dcan be easily adjusted to the 
different values of the parameter p. For this purpose numbers of channel-type functions and 
hyperspherical harmonics have been changed with changing p. It is necessary to emphasize 
that the dependence of the numbers of the trial functions on the parameter p has not been 
exploited in analogous calculations. This dependence gives rise to more flexibility of the 
basis and allows one to avoid numerical instabilities when solving equation 

As a result of the solution of equation (|]) eigenpotentials en{p), Quip), Pi2{p) have 
been obtained. The properties of mesomolecules and transition rates are mostly defined 
by the specific form of the effective potentials £:„(p). The lowest effective potential ei{p) 
describes asymptotically the decay channel H + Hep and is repulsive at all p values. The 
next effective potential £2(p) describes asymptotically the channel He + Hp. As it was 
already mentioned this potential has an attractive part and supports the resonant state we 
are interested. Quip) and Puip) give rise coupling of channels. In case of small coupling 
energy levels of '^Hedp and ^Hedp will be found as eigenvalues of the equation: 

- ^ - ^2(P) - P22(P) + E]u2ip) = (18) 

for zero boundary conditions 
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M2(0) = M2(00) = 0. (19) 

Analogously the continuum wave-function has been found in the one-level approximation as 
a solution of the equation: 

- ^ - ^i(p) - Pnip) + E]uM = (20) 
for the following boundary and asymptotic conditions: 

nifc(O) = 0, 

(21) 

uik{p) — ^ sin{kp + 6), 

p— >cx> 

where k = [E — £1(00))^/^ and phase 5 is of no interest for our purposes. The radiationless 
decay rate is given by 

A = ^|M,^^s-\ (22) 

where the matrix element of the channel coupling operator is 

00 (id 
Mfc = 1 dpuMlQMj-^ + JpQM - Pi2{p)]u2{p). (23) 

III. NUMERICAL RESULTS 

The following values of the masses were used in calculations: = 206.769me, = 
3670.481me, m4He = 7294.295me, ms/fe = 5495.881me. Equation (H) has been solved for 
a number of p values in the interval < p < 45. Variations of the upper bound of this 
interval do not change final results. Expressions (p!3|)-(p!5|) have been used to calculate 
£n{p),Qni{p), Pni{p) for thcsc p valucs. The set of trial functions (|16|) has been adjusted in 
the following way: numbers of channel-type functions A^i and hyperspherical harmonics N2 
were chosen as: 

N-^ = 2, A^2 = 105, p < 5; 
Ni = 6, N2 = 91, 5 < p < 7; 

(24) 

Ni = 6, N2 = 105, 7 < p < 15; 
A^i = 12, N2 = 78, 15 < p. 



The relative accuracy of two lowest eigenpotentials ei(p), £2(p) calculated in the above men- 
tioned interval of p can be estimated as 10""^. Mesomolecular binding energies and 
radiationless decay rates for angular momentum values L = 0, 1, 2 have been calculated 
as described in the previous section. 

The integrand in (^) contains the rapidly oscillating function uik{p), the sharp functions 
U2{p), Quip), Puip) and their derivatives. In consequence of these facts, special care has 
been taken of the evaluation of this integral. For this purpose, U2{p),Qi2{p) and Puip) 
were expressed as a product of sharp functions given in the analytical form and of smooth 
functions given numerically. A few per cent variation of decay rates was found when using 
different ways for analytical representation of sharp functions. 

The calculated values of the binding energies Ebl = — and decay rates for the 
^'"^Hedp systems are presented in Table 1 in comparison with the results of other authors. 

Table 1 
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1 

L .00 
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1 67 
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6Q Q6 


70 6 
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47 QD 


46 7'^ 

4:U. / U 


48 9 
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46 31 
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Q 6 


7 1 1 
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An 








3 "18 


2 87 


Ai 


5.06 




7.0 


2.77 


3.22 


A2 








1.54 


1.74 



Table 1. Binding energies Ebl (eV^) and decay rates (lO^^s ^) of the systems ^'^Hedfi 
calculated in Ref. [0, [0, |TT]] and in the present paper. 



IV. DISCUSSION 



From Table 1 it is clear that binding energies for a given L are close to each other in all 
calculations. One can see that energies of the present paper are higher in comparison with 



calculations |0] and The method of this work gives an upper bound of eigenenergy if 
the coupling of channels is omitted. One can conclude that this fact supports the validity 
of the one-level approximation in our approach. 

The comparison with the results obtained in the framework of the Born-Oppenheimer 
approximation ( !§, |P) cannot be done straightforwardly due to the following reasons. 
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First, mass values and thresholds are introduced in these calculations ad hoc and do not 
coincide with the physical ones. The importance of these procedures for the calculation of 
the decay rate is not clear. Unlike the eigenenergy problem, the calculation of the decay 
rate is very sensitive to the fine details of wave-functions, as is clear from expression (p3[). 
The quasiclassical approximation used in the calculation of the decay rate in the paper ( [Q) 
can be an origin of an additional uncertainty. 

Qualitatively, all calculations support the strong isotopic dependence of the decay rates 
observed in experiment ( [H) Nevertheless, the calculated values are quite different and 
consistency of theoretical results should be reached. 

It is accepted that the formation of Hed/j, molecules takes place in the state with L = 
1. In this connection, for comparison with experiment, the most important is the ratio 
A^/ (A^ + Ai), where A^ is the radiative decay rate from the molecular state L = 1. Using A^ 
from the paper [|lT| and the present values of Ai, one comes to the ratio A^/ (A^ + Ai) = 0.585 
for ^Hedfi and A^/(A^ + Ai) = 0.325 for the ^Hedfi systems. Other processes, which may be 
important in the experiment, are coUisional transitions to the muonic molecular states with 
angular momentum L ^ 1. 

Finally, one would like to emphasize the necessity of the systematic study in the frame- 
work of the same approach of the processes involved in the formation, rearrangement and 
decay of systems under consideration. 
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